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Hydrostatic Journal Bearings (HJBs) are reliable and resilient fluid film rotor support elements 
ideal to replace roller bearings in cryogenic turbomachinery. HJBs will be used for primary space-power 
applications due to their long lifetime, low friction and wear, large load capacity, large direct stiffness, 
and damping force coefficients. An analysis for the performance characteristics of turbulent flow, orifice 
compensated, spherical hydrostatic journal bearings (HJBs) is presented. Spherical bearings allow 
tolerance for shaft misalignment without force performance degradation and have also the ability to 
support axial loads. The spherical HJB combines these advantages to provide a bearing design which 
could be used efficiently on high performance turbomachinery. 

The motion of a barotropic liquid on the thin film bearing lands is described by bulk-flow mass 
and momentum equations. These equations are solved numerically using an efficient CFD method. 
Numerical predictions of load capacity and force coefficients for a 6 recess, spherical HJB in a LCX 
environment are presented. Fluid film axial forces and force coefficients of a magnitude about 20% of 
the radial load capacity are predicted for the case analyzed. Fluid inertia effects, advective and 
centrifugal, are found to affect greatly the static and dynamic force performance of the bearing studied. 


NOMENCLATURE 

A„ Cd 7T d 0 */4. Equivalent orifice area | nr], 

A r 1 R. 0 r . Recess area |nr|. 

C o (S),C.,C 0 Radial clearance function, characteristic clearance [m], CJC, 


C a4t Force damping coefficients due to displacements [Ns/m], ( a,(i = X,Y,Z) 

D. 2 R„ Bearing diameter at midplane (S = 0) [m] 

d„ C d Orifice diameter [m|. Orifice discharge coefficient 

fj.B a M [1 + (c M r J B /H + b M /R J B ) eM ] ; a M = 0.001375, b M = 500,000 

e M = 1/3.00, c M = 10,000, 
r = surface roughness 

Turbulent flow friction factors at journal and bearing surfaces. 

Fx,F Y ,F z Film forces along {X,Y,Z} axes [N] 

f'x?fV?fz cosy, cosy , -siny 

h c H 0 /c. = C 0 (s) + £ x„fxhx + fYolyhy + fzofzhz Dimensionless zeroth-order film thickness 

h, Ae„ f„ h„ . First-order film thickness function. 

bx;by,h z cos0, sin0, +1. Circumferential film thickness components. 

H r Recess depth [m] 

Force stiffness coefficients due to displacements TN/m], (a,0 = X,Y,Z) 

D,L R ,L, Bearing axial length = 1^+ L, [m], Right and left axial side lengths measured from 
recess center 

1 Recess path length [m] 


M U . Circumferential velocity Mach number 

M, yJ} Force inertia coefficients due to displacements [Ns 2 /m], (a,/3 = X,Y,Z) 
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V 

e 

e r 

*.=*« 

«j,*b 

*r 

P,P* 


£iL*£»R 

n,oj 

T 

r, 


Number of recesses on bearing 

Fluid pressure, supply and recess pressures [N/m 2 ] 

Pressures just before and after recess edge [N/m 2 ] 

Discharge pressures on left and right sides of bearing [N/m 2 ] 

Min{P L ,PR}- Characteristic discharge pressure [N/m 2 ] 

(P-P.)/(P.-P»). Dimensionless fluid film pressure 
Dimensionless dynamic (first-order) pressures 
Orifice mass flow rate [kg/s]. 

Bearing radius, Characteristic bearing radius [m], R(s)/R. 

(p Ci C R /p),. Nominal circumferential flow Reynolds number 
(p U C /p),. (C/R).. Modified pressure flow Reynolds number 
(p u C 7 lp)». Nominal Squeeze film Reynolds number 
ip U, C Ip),. Path flow Reynolds number 

(P/P)# )j[U 6 -Q K ) 2 + U‘ (p/p ) H {[ U;+Ul) 

Flow Reynolds numbers relative to journal and bearing surfaces 
Path coordinate on plane of bearing surface [m], S/R. 

Bearing path lengths on right and left sides of bearing [m] 

Fluid mean operating temperature [°K] 

C 2 .(P,-P.)/( m .R.), Characteristic pressure flow speed [m/s] 

(U.,U,)/U.. Dimensionless mean flow velocities in path (s) and circumferential (&) 
directions 

(H r +H)A r +V,. Total recess volume, Volume of orifice supply line [m 3 ] 

Inertial coordinate system 

Z(s)/R,. Dimensionless axial coordinate 

( llp)(dp/dP ). Liquid compressibility coefficient [nr/N] 

( e x» e v ,^z)IC,. Dimensionless journal eccentricities in X, Y, Z directions 

Dimensionless dynamic (perturbed) eccentricities 

Local slope of path coordinate (S) relative to Z axis 

H/(H r +H). Ratio of land film thickness to recess depth 

Circumferential or angular coordinate 

Recess angular length [rad] 

1/2(k,+* b ). Turbulence shear factors in (s,0) flow directions 
f, R J? f B R„. Turbulent shear parameters at journal and bearing surfaces 
(Re h r ) 06SI /7.753. Turbulent shear flow parameter at recess 
Fluid density [kg/m 3 ], characteristic density [kg/m 3 ] 

Fluid viscosity [Ns/m 2 ], characteristic viscosity [Ns/m 2 ] 

Empirical recess-edge entrance loss coefficients in circumferential 
(upstream, downstream) direction 

Empirical recess-edge entrance loss coefficients in path direction (left and right of recess) 
Rotational speed of journal, excitation or whirl frequency [1/s] 

«t. Dimensionless time coordinate 
Recess boundary with outward normal n. 


386 



INTRODUCTION 


Hydrostatic Journal Bearings (HJBs) are the ideal candidates to replace roller bearings as support 
elements in cryogenic turbomachinery. These bearings will be used tor primary space-power applications 
due to their mechanical simplicity, long lifetime, low friction and wear, significant load capacity, and 
large direct stiffness as well as damping force coefficients. HJBs, unlike rolling element bearings, have 
no limit on a DN constraint, and shaft speeds can be allowed to increase to a level more suitable for high 
operating efficiency with reduced overall turbo machinery weight and size. Durability in HJBs is assured 
by the absence of contact between static and moving parts during steady-state operation, while long life 
reduces the frequency of required overhauls. Despite these attractive features, fluid film bearing stability 
considerations due to hydrodynamic and liquid compressibility effects are a primary concern for operation 
at high rotational speeds along with large pressure differentials. The present technological needs call for 
reliable and resilient fluid film bearing designs which provide maximum operating life with optimum - 
controllable rotordynamic characteristics at the lowest cost (Scharrer, 1991, 1992a). 

The analysis of the flow and force response in turbulent flow hydrostatic bearings is of complex nature 
and confined within the realm of classical lubrication theory until recently, see for example Redecliff and 
Vohr (1969), Artiles et al. (1982), Chaomleffel et al. (1986). HJBs for process liquid applications present 
unique flow conditions, requiring for low viscosity liquids large levels of external pressurization to 
provide adequate load capacity and radial stiffness support. Typical pressure drops across a HJB can 
be as large as 30MPa and determine a fully inertial - turbulent fluid flow with significant variation of the 
liquid material properties across the flow region. 

San Andres (1990,1992a) introduced turbulent - inertial bulk-flow models for the analysis of compressible 
liquid (harotropic) HJBs. Extensive numerical predictions have revealed the importance of fluid inertia 
at the film lands and at the recess boundaries of typical high speed HJBs. San Andres (1991a) shows that 
moderate to large journal eccentricities have a pronounced effect on the force coefficients of HJBs with 
large hydrodynamic effects (high rotational speeds). Furthermore, orifice back-flow along with a sudden 
drop on direct stiffness are likely to occur at large eccentricity operation. 

Kurtin et al. (1991), Franchek (1992), and Mosher (1993) present relevant experimental data for the static 
and dynamic force characteristics of water lubricated, turbulent flow, hydrostatic bearings. Experimental 
measurements are routinely performed for hydrostatic bearings of different geometries and at journal 
speeds ranging from 10,200 to 24,600 rpm and pressure supplies from 4 to 7 MPa. These references 
also present extensive comparisons of test results with numerical predictions based on the models of San 
Andres (1990, 1992a). In general the correlation between experimental and theoretical results is very 
good for conventional HJB geometries. It is noted that accurate theoretical results depended greatly on 
the knowledge of the bearing operating clearance, and most importantly, on the orifice discharge 
coefficients. 


Adams et al. (1992) have also presented test results for the rotordynamic force coefficients of a four pad, 
one recess/pad, laminar flow, hydrostatic bearing. The experiments were performed with SAE 30 oil 
and at low rotational speeds (1000 and 2000 rpm) and low pressure supplies (max. 2.6 MPa, 375psig). 
Force stiffnesses and direct damping coefficients seem to be well identified while cross-coupled damping 
and inertia force coefficients show a rather unexpected behavior. 

The threshold speed of instability and the whirl frequency ratio (WFR) define the stability 
characteristics of a simple rotor-bearing system. This instability is of the hydrodynamic type and solely 
due to the effect of journal rotational speed on the flow field. Incompressible liquid hydrostatic bearings 
present a whirl frequency ratio identical to that of plain journal bearings (WFR -0.5). This condition, 
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as also verified experimentally by Mosher (1993), then limits severely the application of HJBs to high 
speed, light weight turbomachinery. HJBs handling highly compressible liquids, such as LH 2 for example, 
are prone to show a self-excited type instability of the pneumatic hammer type and can produce negative 
damping force coefficients for low frequency excitations. Dynamic operation under these conditions will 
then result in a poor stability indicator (WFR) greater than 0.50. This important result, although first 
reported by Redecliff and Vohr (1969), has been largely overlooked until recently. 

Recommended fixes to improve the limited dynamic stability of turbulent flow hydrostatic 
bearings are: 

• Use of large scale roughened bearing surfaces to reduce the cross-coupled stiffness coefficients 
directly promoting hydrodynamic bearing instability (Franchek,1992). 

• Use of end seal restrictions or wear end-rings (Scharrer et al., 1992b) to control bearing leakage, 
increase the damping coefficients, and add a degree of safety for start-up and shut-down transient 
operation. 

• Use of liquid injection opposing journal rotation to reduce the development of the circumferential 
flow velocity and eliminate the cross-coupled stiffness coefficient (Franchek, 1992). 

The development of a leading technology in HJBs calls also for a bearing geometry not only able 
to provide radial load support but also with the capability to handle axial loads accompanied by shaft 
dynamic axial excurtions (Sutton et al., 1991). Spherical bearings offer this advantage along with the 
capability to tolerate large levels of static and dynamic misalignment (from journal and bearing) without 
alteration of the bearing performance. Furthermore, the recent developments in CNC manufacturing 
processes allow to machine spherical surfaces almost as quickly and economically as cylindrical surfaces 
(Craighead, 1992). 

Goenka et al. (1980) and Craighead et al. (1992) have provided analysis of spherical journal bearings 
for laminar/turbulent flow applications. However, in high-speed, turbulent flow applications with process 
liquids of low viscosity, fluid inertia effects need to be accounted. Most notably, centrifugal flow 
acceleration terms are of particular importance for these operating conditions. The present study 
considers the analysis of hemispherical hydrostatic bearings with a barotropic liquid. Turbulent bulk-flow 
equations of motion are derived and solved numerically using an efficient CFD algorithm. Numerical 
predictions for the load capacity (radial and axial) and dynamic force coefficients for a L0 2 HJB are 
presented and discussed in detail. 


ANALYSIS 

Consider, as shown in Figure 1, the flow of a variable properties liquid in the thin film region 
between an inner rotating journal and a stationary bearing. Cryogenic liquids are characterized by low 
viscosities, and thermal (energy transport) effects due to friction heating and kinetic energy variations are 
expected to be of minor importance in the performance of hydrostatic bearings. This assertion is not fully 
justified for especial operating conditions (see for example Yang et al., 1992a). On the other hand, due 
to the large levels of pressure differential required to provide substantial load capacity, the effects of 
pressure on the liquid properties, and ultimately on bearing performance, are thought to be of primary 
importance. 

Figure 2 shows the journal outer surface as a surface of revolution formed by rotating the curve 
R(Z) about the axis Z. The path (S) and circumferential ( 0 ) coordinates are used as independent spatial 
variables. The coordinates {Z.R} defining the journal surface are expressed as parametric functions of 
the path cootdinate S. Trigonometric function of the angle y defining the local slope of the path relative 
to the axis (Z) are given by: 
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( 1 ) 


tan 7 


dR. 

dZ' 


sin 7 = - 


dR, 
dS ’ 


cos 7 = 


dZ 

dS 


where the coordinate relationships for journal/bearing spherical surfaces are: 

7 = R(S) = fl.cos 7 , Z(5) = R t sin 7 

and R.(=D./ 2 ) corresponds to the journal radius at the bearing axial midplane. 


( 2 ) 


At operating conditions, the journal position relative to the bearing housing is described with 
reference to the inertial axes {X,Y,Z} by the journal center displacements (ex(t), eyft), e^t)). Simple 
geometrical relationships determine the film thickness in the flow region to be given by the following 
expression (Goenka et al., 1980, Childs, 1989): 


H(S,6,t ) = C(S) + e x cos 7 cos 0 + ^ cos 7 sin 0 - e z sin 7 ( 3 ) 

In a spherical bearing, journal axis rotations or misalignment provide no film thickness variation on 
the flow region. C(S) in (3) above is a general function describing the radial clearance variation along 
the path coordinate for the journal centered position. 


The equations of motion 

The turbulent bulk-flow equations for a variable properties (barotropic) liquid on the thin film lands 
of the spherical bearing are given as (Childs, 1989): 


Equation of continuity: 




(4) 


Path momentum equation 

h E - ji (ft u) + 
dS H ' 1 dt 


+ 

R 


d(pHU s UJi) d(pHU e UJ rjrj2 dR 

dS— + —He— - pHU • Ts 


(5) 


Circumferential momentum equation: 


-H— = 

Rdd 



d(pHU t ) 1 

dt R 


dipHUJUJQ ^ d(pHU e U $ ) 

ds + de 


+ P HU t U, *!L 
’ * dS 


( 6 ) 


on the region {-S, <S<S R ; O< 0 < 2 tt}; and where, k,=k # =(k 1 +k B ) are the wall shear stress difference 
coefficients taken as local functions of the turbulent friction factors, Reynolds numbers and surface 
conditions, i.e, k J =f J R J , k B =f B R B (Hirs, 1973, San Andres, 1992a). For inertialess-laminar fluid flows 
the equations above reduce to the classical form given by Goenka and Booker (1980) for spherical bearing 
geometries. 


For cryogenic liquids such as LH2, L02, LN2, and LCH4, the fluid properties are calculated from 
the Benedict-Web-Rubin equation of state as given in the standard computer program and data base of 
McCarty(1986). 
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Recess Flow and Pressure equations: 

A mass conservation equation at each bearing recess of area (I. R..0 r ) and depth H r is defined by the 
global balance between the mass flow through the orifice restrictor (Q ro ), the mass flow into the film lands 
and the time rate of change of liquid mass within the recess volume V,. This equation is given as: 

. ey/1: (7) 

for r = 1,2,..., Nrecess 


where /3=(l/p)3p/dP represents the fluid compressibility material coefficient at the recess volume, and 
T r is the closure of the recess volume with the film lands and with normal n along the boundary line. 
Note that the orifice flow equation is valid only for small changes of the liquid density (Hall et al., 
1986). 

The fluid edge pressure at the entrance to the film lands is given by the superposition of viscous 
shear effects on the recess extent and an entrance drop due to fluid inertia. On the circumferential 
direction, the pressure rise (P c ) downstream of the recess orifice is given by (Constantinescu et al., 1987, 
San Andres, 1992a): 


P; = P r - VrK 


*.e 

2H; 


U e (p' t lp r )TI - 


QR(s) 


(1 - M 2 ) 


( 8 ) 


where M is the circumferential flow local Mach number at the orifice discharge. 

The entrance pressures (P e ) to the film lands in the circumferential and axial directions are given 
by: 


p; = p; - y(i+q{i - (p>;)V} Vl 


(9) 


P: = Pr - y (1 

for r = 1,2.,., Nrecess 

The Bernoulli like pressure drop in equations (10) is considered only if the fluid leaves the recess 
towards the film lands. If on the contrary, fluid enters from the film lands into the bearing recess, then 
the edge pressure takes the value of the recess pressure (P r ). This consideration is based on momentum 
conservation for turbulent shear flows in sudden expansions and also on the fundamental measurements 
of Chaomleffel et al.(1986). The inertial pressure drop given above does not account for centrifugal flow 
effects in the spherical bearing geometry since the change in the bearing radial coordinate from recess 
edge to film lands is small. 


*,){i - (p;i P ;) 2 r, 2 }u; 


( 10 ) 
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Boundary Conditions: 


Due to periodicity, the pressure and velocities are continuous and single-valued in the circumferential 
direction 0; i.e.. 


P, U „ U ,( S , Q,t ) = P, U s , U ,( S,e + 2tt,0 (11) 

At the bearing side discharge planes, the fluid pressure is equal to specified values of discharge or 
sump pressures, i.e.: 

at the right plane , Z = +L K P(+S H ,6) = P H (d) 

at the left plane, Z = -L, P(-S L ,0 = P H (0) (12) 

In general the discharge pressures are uniform and constant. However, in some cryogenic turbopump 
applications the bearing may be located close to the pump-impeller discharge. In this case, the sump 
pressures are non-uniform though rotationally symmetric and expressed by a Fourier series. The boundary 
conditions described are valid for fluid flows well below sonic conditions. 

Perturbation Analysis 

Consider the journal center to describe small amplitude harmonic motions about an equilibrium static 
position. That is, let the journal center displacements be given as 


e xW = e Xo * Ae x e'»', eft) = e Ko + Ae y e eft) = + Ae z e™'; i =fT (13) 

where « denotes the frequency of the whirl motion. The magnitudes of the dynamic perturbations in 
journal displacements, l{Ae x ,Ae v ,Aez}/CJ are very small (i.e. < < < 1). Then, the film thickness is 
given by the superposition of steady-state (h u ) and dynamic (h,) components given by the real part of the 
expression: 


h = h 0 + h l e h " (14a) 

where h {t = C^(s) + {c Xo cost? + e Yo sin0} cosy - c. /M sin-y (14b) 

h t = K = Af x/x K + A tyfy K + A f 2 /z h z (14c) 

with ffs) = ffs) = cosy; f 2 = -siny, and h x = cos0; h r = sin0, h z = 1 (15) 


are the film thickness perturbed functions along the path and circumferential coordinates, respectively. 
These functions greatly facil itate the comprehension of the perturbed flow field equations and the resulting 
rotordynamic force coefficients given latter. 

The flow field variables (U„U 0 ,P), as well as the fluid properties (p,p ) and the shear parameters 
(ke*k») ar £ also formulated as the superposition of zeroth-order and first-order complex fields describing 
the static equilibrium condition and the perturbed dynamic motion, respectively. In general, these fields 
are expressed as: 
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* = + ^"{ Ac A + Ac A + Ae z ^z} = + ‘''Ac.*. , 


« = XJ,Z 


(16) 


San Andres (1990, 1992a) and Yang (1992b) discuss the procedure for the numerical solution of the 
non-linear flow equations. The differential equations of motion are integrated on staggered control 
volumes for each primitive variable. Program computing time is relatively small since the code uses 
accurate approximate analytical solutions to initiate the computational procedure and accelerate 
convergence to a solution defined by the operating parameters and bearing geometry. 

Fluid Film Forces and Dynamic Force Coefficients 

Fluid film forces are calculated by integration of the zeroth-order pressure field on the journal surface, 


>*■ 

= (p,-p.)R: \‘‘ 

[ 0< * ? ' p • 

f x • COS0 

/• sin0 

1 

J -h 


1 

1 


(17) 


The perturbation analysis allows the dynamic force coefficients due to journal center displacements 
to be obtained from the general expression for dynamic forces given as: 




~ K K K 

iX xx XY ^xz 
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'**x~ 

af y 

= - 
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c c c 
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K ZX K 7J K 7Z_ 
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c c c 

w zx ^ / r ^ zz 




M yjx M zy 


My._ 


(18) 


where the force coefficients due to journal center displacements are given by: 


- <•>’ * ’ C* 


, . (r.-r.yt: 


r s. r 9 . * 7* 

l: Ptf.-Kr-ti 


cid 


(19) 


= X,Y,Z 


RESULTS AND DISCUSSION 

Experimental results for the steady state and dynamic force response characteristics of turbulent flow 
HJBs for process liquid applications are given by Kurtin et.al (1991), Adams et.al (1992), Franchek 
(1992), and Mosher (1993). These studies are relevant to the investigation of cylindrical bearing 
geometries with large pressure drops and high rotational speeds similar to those found in high 
performance turbomachinery components. Correlation of test measurements with predictions based on 
the present flow model are very favorable for smooth surface HJBs (Franchek. 1992, Yang. 1992b 
Mosher, 1993). 
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Spherical hydrostatic bearings may provide a unique alternative for radial load support in cryogenic 
liquid turbopumps (Sutton et.al, 1991) since they also offer the distinct advantages of tolerance to journal 
misalignment and ability to withstand shaft axial motions by providing axial thrust. In the following, the 
static and dynamic performance characteristics of a spherical HJB geometry handling LO ; are presented 
and discussed. Normalization of results and extensive parametric studies accounting for variations in the 
bearing geometrical and operating parameters would be impractical due to the complex nature of the flow 
field and force response in turbulent hydrostatic bearings. It suffices to say that the example presented 
corresponds to a bearing element designed for optimal radial support at the rated operating conditions. 

Table 1 shows the geometry and operating conditions of a 6 recess, LO ; hydrostatic bearing with 
a fixed pressure drop across the bearing (2,000 psi) and a rotational speed equal to 22.5 Kcpm. The 
spherical bearing diameter (D.) and axial length (L) are equal to 91.036 mm and 64.37 mm, respectively. 
The exit diameter of this bearing is equal to 64.37 mm and the arc described by the spherical path 
between the bearing middle and discharge planes is equal to 7 *= 45 °. At the rated conditions, a recess 
pressure ration (P r ) equal to 0.55 provides maximum direct stiffness coefficients and requires orifices of 
diameter equal to 2.37 mm. Table 1 also includes values for the empirical recess-edge non - isentropic 
loss parameters (£) and orifice discharge coefficients (C d ) used in the analysis. The values chosen are 
representative from those used in the extensive experimental - theoretical studies of Kurtin et al. (1991) 
and Franchek (1992). The regime of operations of the bearings is fully turbulent with circumferencial 
(Rec) and axial flow (R ai ) Reynolds numbers equal to 56,900 and 35,857, respectively. A comprehensive 
study and comparison between the performance characteristics of equivalent cylindrical and spherical 
HJBs can be found in the work of San Andres (1992b). 

Numerical predictions are presented for the static and dynamic force characteristics of the bearing 
for increasing values of the axial journal eccentricity (e z ) while the journal center is displaced radially 
towards the middle of the bottom recess, i.e. e x varies and e Y =0. From equation (3), the maximum 
axial journal displacement is equal to & z = ( 1 - e x cos 7 ’ )/sin 7 *, where 7 * corresponds to the spherical 
angle at the bearing discharge, that is 45° for the example presented. 

Figure 3 shows the radial load of the bearing as the static eccentricity ( e^) increases and for values 
of axial journal displacement (e^ equal to 0.0 and 0.60, respectively. The results show the load to 
increase linearly with the journal lateral displacement denoting a bearing with uniform stiffness 
characteristics for most eccentricities. The effect of the axial journal motion is relatively small on the 
total bearing load. Figure 4 shows the restoring force (-F z ) as the axial journal is displaced towards the 
bearing shell and for increasing values of the static lateral eccentricity (c x ). The axial force appears to 
be linear with displacement and increases with the radial journal displacement. Note that the thrust load 
(F z ) is about 1/5 of the radial load, and shows the spherical bearing to have a limited axial load capacity 
in comparison with its radial load support. 

Figures 5 and 6 show the direct (K xx ,Kyy) and cross-coupled (K^, K yx ) radial force stiffness 
coefficients, and Figure 7 presents the direct damping coefficients (C xx ,Cyy) for increasing values of the 
journal eccentricity e x . The figures show the axial journal center displacements (e z ) not to affect these 
coefficients except at large lateral eccentricities (c x ). Note that the radial force coefficients are relatively 
constant for radial eccentricities as large as 50% of the bearing clearance and show clearly the major 
benefit of a hydrostatic bearing. A lucid discussion on the effect of these radial coefficients on the 
rotordynamic lateral force response can be found elsewhere (San Andres, 1991a). 

Figures 8 and 9 show the direct axial stiffness (K zz ) and direct damping (C zz ) coefficients for 
dynamic journal axial motions as the lateral eccentricity (e x ) increases. These coefficients increase with 
the journal axial position and show a significant rise at moderately large radial eccentricities. At the 
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concentric position (c x = e Y = 0), the axial stiffness is about 1/5 of the radial stiffness (K xx ), and the 
axial damping C TJ is approximately 11% of the direct lateral damping (C xx ). Thus, the axial force 
coefficients are relatively small when compared to the radial force coefficients. On the other hand, the 
magnitudes of the axial force coefficients are still significant in terms of their ability to sustain limited 
dynamic axial load conditions. For example, from Figure 4 it is inferred that the spherical HJB can 
tolerate safely an axial load as large as 5,000 N (1,125 lbs). 

Cross-coupled axial force coefficients (K 7X ,K ZY ,C ZX ,C ZY ) due to journal lateral motions, and radial 
force coefficients (K XZ ,K YZ ,C XZ> C YZ ) due to journal axial motions are not reproduced here for brevity. 
The radial force coefficients are very small except for large journal center displacements and provide a 
measure in uncoupling between axial dynamic motions and lateral dynamic force response. The axial 
force stiffness (K ZX ,K ZY ) are approximately 1/2 of the direct stiffness (K zz ) at e z =0.80, e x =0, and 
decrease rapidly as the lateral eccentricity (e x ) increases. 

For completeness in the analysis, calculations were performed to determine whether fluid inertia 
effects, both advective and centrifugal, are of importance on the static and dynamic force performance 
characteristics of the spherical hydrostatic bearing. Table 2 presents a summary of the results for 
increasing values of the journal rotational speed at the bearing concentric position ( e x = c y = c z =0). 
The orifice diameter and loss coefficients are identical for the simulations. A comparison of results 
shows large differences on the recess pressure ratio (P r ) and the force coefficients. Fluid inertia acts as 
an additional flow resistance and then determines larger recess pressures with a reduced flow rate. The 
most notable effect is related to the reduced magnitude of the direct radial stiffnesses (K xx = Kyy). At 
the rated operating point, 22.5 Kcpm, the values of direct stiffness are equal to 240.9 and 308.5 MN/m 
for the bearing with and without fluid inertia effects. This corresponds to a net reduction in hydrostatic 
load capacity of 22% and it is a direct consequence of the increased flow resistance due to centrifugal 
effects on the curved flow path and also due to the larger recess pressures. It is evident from the results 
presented that fluid inertia effects need to be included in the analysis of high speed, turbulent flow HJBs. 

CONCLUSIONS 

An analysis for the performance characteristics of turbulent flow, orifice compensated, spherical 
hydrostatic journal bearing (HJBs) is presented. Hydrostatic bearings offer a substantial radial load 
capacity and can be used with process liquids of low viscosity if large pressure differentials across the 
bearing are available. On the other hand, the spherical bearing geometry allows tolerance for shaft 
misalignment without force performance degradation and it also has the ability to support thrust loads. 
The spherical HJB combines these advantages to provide a bearing design which could be used efficiently 
on high performance turbomachinery. 

The motion of a barotropic liquid on the thin film bearing lands is described by bulk-flow mass and 
momentum equations. Zeroth-order equations describe the fluid flow field for a journal equilibrium 
position, while first-order linear equations govern the fluid flow for small amplitude journal center radial 
and axial motions. Solution to the zeroth-order flow field equations provides the bearing flow rate, 
radial/axial film forces and drag torque. Solutions to the first-order equations determine the rotordynamic 
force coefficients due to journal lateral and axial motions. Numerical predictions of load capacity and 
force coefficients for a 6 recess, spherical HJB in a LO : environment are presented for increasing values 
of the journal center radial and axial displacements. The results show that axial journal motions do not 
alter significantly the radial load capacity of the bearing. On the other hand, the spherical bearing 
geometry provides fluid film axial forces of a magnitude about 20% of the radial load capacity for the 
example analyzed. Fluid inertia effects, advective and centrifugal, are found to affect greatly the static 
and dynamic force performance of the bearing studied. 
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Table 1 . Geometry and Operating Conditions of 6 recess. LO 2 
Spherical Hydrostatic Journal Bearing. 

Dimensions : Numner ot recesses Nrec - 6. 

Clearance, c-101 ,6umi0.004 in), recess oeoin. Hr-381 urn, 0.01 5in): H«c=3.75 

-pOencal: Diameter. D*-91 036mm (3.58mi: Path Lengin. S -71 50mm (3.315in) 
Dr- 64.372mm (2. 5534ml; 

Axial lengtn. L- 64.372mm (2.5534m): LD‘ - 0.707. R*/c-448.0 
Spherical angle yil « 45 deg. 

Recess: length /=31.74mm( 1.25 in), circumierential lengtn 0r=3O.O oeg 
journal and bearing sunace conditions: smooth 

Orifice Ca-0.89. ...ameter do*2. 377mm :z r concentric recess cresure ratio pr=0.55, 

Recess edge (non-isentrooici coerticients c^-O.O; Co.0.50; ',-=00 

9pe rating Parameters: ., rotational speed: 2.356 rad/s (22.5 Kcpm) 

Pressure supply. Ps-15.16 MPa (2.200 psiai 
exit. Pa- 1.38 MPa ( 200 psiaj 

F]ui^ LO 2 at 90K r!52 R). Us=0.22454E:V3 Pa-s. ps- 1 .172 Kg/m3 
,Ua*0.19813E-3 Pa-s. pa-1.144 Kg/m3 
(3 -0.1696-8 1/Pa-1„ (85.52 kpsi).. 

TYP Reynoids numbers: 

Rectrc»0si)R*C/Us=o6.9OO ; 

Reu.aj-rrv(2TiDR Uii= 35.857. Mass How rate concentnc-2.3735 kgys 




recess orifice 


LV'2 




Soherical ; ^ 0 

HJ5 ^ " 



Table 2 

Spherical hydrostatic bearing performance characteristics 
Effect of Fluid Inertia 

(a) Model vith fluid inertia effects, d *2. 37mm 

o 

Speed Pr-Pa Mass flow Torque ^ C^-C^ C, 


(Kcpm) 


(kg/s) (N»a) 


(MW/m) 


0.00 0.4700 3.1314 0.000 179.32 0.00 6.82 151.90 0.00 19.30 

12.50 0.4980 3.0160 2.833 229.20 112.30 6.74 153.60 16.86 19.50 

22.50 0.5000 2.8735 5.980 240.90 198.90 6.43 165.30 25.50 19.70 


(b) Model without fluid inertia effects, d *2. 37mm 
o 

Speed Pr-Pa Mass flow Torque ^ C^-C^ C, 


(Rcpa) a ~ ra (kg/s) (N'b) 

0.00 0.4157 3.2763 0.000 257.60 

12.50 0.4302 3.2290 2.996 280.40 

22.50 0.4632 3.1259 6.086 308.50 


(MM/e) (MM-s/a) 

0.00 5.07 92.38 0.00 10.40 
39.40 4.88 101.30 2.10 10.30 
64.09 4.71 110.50 2.24 10.30 
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Figure 2. Coordinate Relationships In a Spherical Bearing 
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FORCE -Fz (KN) 



RADIAL STATIC ECCENTRICITY (e,) 

Figure 3. Radial Load vs. Journal radial eccentricity (Ex) for Journal axial 
displacements Ex=0.0 and 0.60 



AXIAL JOURNAL DISPLACEMENT (e*) 

Figure 4. Axial force (-Fz) vs. Journal axial eccentricity (Ez) for Increasing 
Journal radial displacements Ex 
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STIFFNESS — Kyx AND Kxy (MN/m) 



Figure 5. Lateral force stiffness coefficients (Kxx, Kyy) vs. Journal radial 
eccentricity (Ex) for Journal axial displacements Ez=0.0 and 0.60 



Figure 6. Lateral force stiffness coefficients (-Kxy, Kyx) vs. Journal radial 
eccentricity (Ex) for Journal axial displacements Ei=0.0 and 0.60 
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DAMPING Cxx AND Cyy (KN-s/m) 


Cxx (e x =0.0) 
Cxx (e x =0.6) 
Cyy (e z =0.0) 
Cyy (e z =0.6) 


/ 


"~o'o • o'. o!2 0:3 O.H 0.0 06 0.7 08 < 

RADIAL STATIC ECCENTRICITY (e,) 

Figure 7 Lateral force damping coefficients (Cxx, Cyy) vs. 

9 eccardrW M lor Joun»l «W dlsplscanrsnts t~0.0 and 0.60 


RADIAL STATIC ECCENTRICITY (e K ) 

FloureS. AM loros stillness coedlclenttKzz) vs. loumsIradM eccentric!,, 
fc x ) f 0 r increasing Journal axial displacements Ex. 


dampin 



RADIAL STATIC ECCENTRICITY (e,) 


Fl 8 u re9 . AM torce damping eoemc^nttC^, vs. /oumalraM eccnMcHy 
(£x) for Increasing Journal axial displacements e* 
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